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ABSTRACT 


For acoustic frequencies of interest in the development of undersea wireless wide- 
area networks, wind noise is the dominant eontribution to channel noise. The large 
dynamic range of wind noise forees the network designer to eonsider wide variations in 
link margin and effeetive node-to-node range. Previous reported eorrelations of aeoustie 
eommunieation performanee and wind speed led to multiple hypotheses explaining 
degradation, ineluding noise variability as well as wind-driven sea-surfaee effeets sueh as 
roughness, entrained bubbles, mixing, and stratifieation. This thesis examines operations 
of a 40-node wide-area aeoustie network in varying noise eonditions. This data set is 
unique in that the environment is strongly downward refraeting. Based on theoretical 
models, wind noise is assumed to dominate all other noise sourees. Considerations 
inelude dependenee of link margin and range on wind speed. Although the experiment 
was not designed to examine the eorrelation between transmission range and wind speed, 
a weak eorrelation was observed in the limited data set available. 
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I. 


INTRODUCTION 


An increasing interest for the United States Navy is the use of underwater 
eommunieations involving High Frequency (HF - greater than 1 kHz) aeousties. This 
teehnology enables the ereation of underwater wireless networks. The performanee of 
sueh a network is dependent on loeal environment and node geometry. But even for a 
given environment and geometry, eommunieation performanee varies with loeal wind 
noise. The working hypothesis is that eommunieations range degrades with inereasing 
wind noise. This thesis attempts to quantify that relationship through the development of 
a link budget model ealibrated with experimentally obtained eommunieations 
performanee data. 

In the next chapter, we develop a link budget model as a form of the sonar 
equation. A link budget is a method used for estimating the signal to noise ratio (SNR) 
for given parameters sueh as transmit power, antenna gain, propagation loss, and 
interferenee. A theoretieal model shows how the link budget changes with varying wind 
values. 

In the third ehapter, we deseribe a relevant experiment. The network design and 
eommunieation are explained. Information speeifie to the network is discussed, sueh as 
the souree level and frequeneies used. 

In the fourth ehapter, we eonsider the environment and the geometry of the 
experiment. 

The fifth chapter graphieally examines how, during the experiment, the wind 
affeeted the performanee of the network. 

In the sixth chapter, we take what was observed in the experiment and re-evaluate 
the link budget model. This provides us a better understanding of how wind affeets 
network performanee. 

Finally, the seventh ehapter eoneludes this thesis with a summary and 
reeommendations for future experimentation. 
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II. LINK BUDGET MODEL 


A. SONAR EQUATION 

In a standard writing of the sonar equation the left hand side is SE, for signal 
exeess. In developing the sonar equation as a link budget model (rather than deteetion 
model), we seek to represent the system’s ability to maintain a eommunieations link in a 
given environment. A eomplete link budget yields the link margin based on signal 
proeessing and other system eonsiderations. In this thesis we limit our development to 
SNR, expressed as: 

SNR = PSL -TL- NSL + (2.1) 

where 

SNR is Signal to Noise Ratio at the reeeiver 

PSL is Pressure Speetrum Level of transmitting platform 

TL is Transmission Loss of the medium 

NSL is Ambient Noise Speetrum Level of the environment 

D1 is Direetivity Index of the reeeiver and transmitter. 

Eaeh term is in deeibels (dB). In this ease, the DI terms are zero, beeause the transmitters 
and reeeivers are omni direetional. Beeause of system and ehannel variation aeross the 
operating speetral band, the remaining terms are a function of acoustic frequency,/. 

1. Pressure Spectrum Level 

The Pressure Spectrum Level {PSL) represents the density of the transmitted 
energy across the spectral band of interest. Lor a transmitted signal composed of tonals, 

P5L = 5L-101ogio(A) (2.2) 

In the calculation, SL is the source level of the wideband transmitted signal, and N 
is the number of tones instantaneously activated. Lor this thesis, we consider M-ary 
Lrequency Shift Keying (MLSK) modulation with 128 tones and l-of-4 keying, or A=32. 
Lor a Source Level of 179 dB re IpPa @ 1 meter, the PSL is 164 dB. 
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2. Transmission Loss 

Transmission Loss {TL) reduces the sound amplitude as the sound propagates 
from the source to a receiver. The reference intensity is measured from a point 1 meter 
from the center of the source. In mathematical terms, TL is expressed asL 

(I ^ 

rL = 101og,o (2.3) 

v-'r y 

In this case, /r is the intensity at the receiver, and /o is the intensity I meter from the 
acoustic center of the source. Transmission Loss is, in general, due to both geometric 
spreading and attenuation, which includes loss mechanisms at the surface as well as at the 
bottom of the sea. 

Transmission Loss is dependent on the local environment, acting on sound energy 
as it propagates along the path from the source to the receiver. TL from geometrical 
spreading is often assumed to correspond to a combination of spherical and cylindrical 
spreading, and is calculated to as: 

TLc^om =101ogio(r,) + 101ogio(r) (2.4) 

Where r, represents the transition range from spherical to cylindrical spreading, 
and r is range, in units of meters. 

In Figure I, the change in the slope is a result of the assumed transition from 
spherical to cylindrical spreading. In this case, the transition range is assumed to equal 
the bottom depth of 220m. 


1 Urick, Robert, (1983). Principles of Underwater Sound Edition (McGraw-Hill, Inc., 1983) 
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Transmission Loss due to spreading 



Figure 1. Transmission Loss vs. Range, where range is logarithmic; the transition 

range is set to 220 m. 

Geometrical spreading is not the only component of TL. As sound propagates, 
there is attenuation of the signal due to absorption, scattering, and leakage. There are 
three effects that contribute to absorption: shear viscosity, ionic relaxation, and relaxation 
time. TL due to absorption varies linearly with range, rather than logarithmically, and 
can be calculated from2: 

TK..en=(^rx\0-^ (2.5) 

^^ 0.11/ 44/ +3 Oxio-4(y2)+33xlO-^ (2.6) 

Where r is range in meters, a is the attenuation coefficient in dB/km, and/is frequency in 
kHz. This adds linearly with the geometric transmission loss, to give a total TL term that 


2 Urick, Robert (1983). Principles of Underwater Sound 3"^. Edition (McGraw-Hill, Inc., 1983). 
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is dependent upon range and frequeney. The effeets of seattering and leakage are 
assumed to be negligible eompared to other loss meehanisms. 

B. AMBIENT NOISE 

Ambient noise is the baekground noise of the oeean due to either natural or man¬ 
made eauses; various sourees eontribute noise in distinet frequeney bands, so that the 
term varies with frequeney. The Noise Speetrum Level (NSL) is a measure of the speetral 
density of the ambient noise in dB. 

1. Wind 

Wind speed is thought to affeet ambient noise in several ways. First, by simply 
blowing aeross the surfaee of the water, the wind ereates flow noise. Also ereated by the 
wind are surfaee waves. When these waves gain large enough amplitude, they make 
noise by breaking on the surfaee, ereating bubbles of air under the surfaee, whieh ean 
oseillate and burst.3 Wind noise is estimated by4: 

=50 + 7.5V^ + 201og,o(/)-401og,o(/+0.4) (2.7) 

where w is wind speed in m/s, and/is frequeney, in kHz. 

2. Precipitation 

Rain falling on the surfaee of the oeean ereates sounds by pereussion, and the 
oseillation and eollapse of bubbles produeed, just as rain falling on a roof during a storm 
makes noise inside a house. Although there were time periods of our experiment where 
preeipitation oeeurred, estimates of times and rates were not available, and the effeets of 
preeipitation are not ineluded in our study. 

3. Thermal Noise 

At high frequencies the noise spectrum is dominated by noise arising from the 
random motion of water molecules. The contribution to the noise spectrum from thermal 
noise, NSLthermai can be calculated firomS; 


3 Urick, Robert (1979) Sound Propagation in the Sea (Defense Advanced Research Projects Agency 
1979). 

4 Wenz, G.M. (1962). “Acoustic Ambient Noise in the Ocean: Spectra and Sources” Journal of the 
Acoustical Society of America, 34,1. 

3 Urick, Robert (1983). Principles of Underwater Sound 5/ Edition (McGraw-Hill, Inc., 1983). 
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"Sin™, =-15 + 201og,„(/) (2.8) 

where/is frequeney in kHz. 

4. Shipping 

Shipping traffie on the open oeean ereates a substantial amount of noise 
underwater, from eomponents sueh as propellers, engines, and flow. This noise ean be 
deteeted over very large distanees. The eontribution from eaeh individual ship may not 
have a large affeet on Noise Level; however, when it is eonsidered that there are 
thousands of tons of shipping erossing the Paeifie Oeean at any time, this is a signifieant 
souree of noise. Aeeuraey of estimates of shipping noise are limited by the faet that 
shipping density and ship meehanies have ehanged in the deeades sinee original estimates 
were made, and vary loeally as well. However, it tends to be below the frequeneies of 
interest for eommunieation, making it less important for our example. Uriek gives the 
following expression for the ambient noise due to shipping. In this expression D is the 
density of shipping traffie, from 0 to 1; 0 being no traffie, and 1 being maximum traffie^. 

NSL^„„ =40 + 20(£)-0.5) + 26log,„(/)-60log,o(/+0.03) (2.9) 

where/is frequeney in kHz. 

5. Biologies 

The animals that live under the surfaee also ereate noise. The snapping of shrimp, 
ealls of whales, and whistles of dolphins are just some of the sounds that are 
eommonplaee, all eontributing to the ambient noise levels. For this analysis, we negleet 
biologieal noise, as it tends to be outside the frequeney range of interest to this study. 

C. AMBIENT NOISE VS. WIND 

All of ambient noise eomponents for whieh speetral level estimates were available 
are eombined into the total NSL and plotted in Figure 2. Clearly, in the frequeneies of 
interest to our study wind noise should dominate the total noise. 


6 Urick, Robert (1986) Ambient Noise in the Sea (Peninsula Publishing Los Altos, CA 1986). 
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Figure 2. Noise Speetrum Level eurves using theoretieal equations; in the HF band, 
wind noise (wind speed 5 m/s, and shipping is medium) is dominant. 

Figure 3 deseribes the variability of the ambient Noise Speetrum Level for 

ehanging wind speeds. Notiee that for HF signals, there is a large dynamie range of more 

than 30 dB. This further indieates the importanee of wind noise in ealeulating ambient 

Noise Speetrum Levels and SNR. 
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Figure 3. Dependence of ambient Noise Spectral Level upon wind speed; note the 
variance in the Fligh Frequency band (greater than 1 kFlz) (dB re 1 pPa). 

D. SNR VS. WIND 

Let us now examine the influence of wind on SNR. Revisiting the sonar equation 

gives: 

SNRa,.„„,,=PSL-TL-NSL, ( 2 . 10 ) 

where the term (-TL-NSL) represents the environmental factor which needs to be 
overcome by the source PSL in order to obtain the SNR necessary to achieve desired 
communication rates. For this study Transmission Loss was assumed to be given by 
geometrical spreading (equation 2.4) with the transition range equal to 220m. 

Attenuation due to absorption was calculated using equation 2.5 and 2.6. The shipping 
level was assumed to be medium, and the sea state is calm. Only the wind speed is 
varied, from 0 to 25 m/s in 5 m/s increments. Figure 4 shows this environmental factor 
for those six different wind conditions. Notice that although the largest impact of 
NSLwind is in the highest frequencies, the largest amount of variability in the impact of the 
wind noise is around 20 kHz, presumably because of the competing factors of wind noise 
and absorption. Wind noise decreases with frequency while absorption increases. This 


9 

























































effect is clearly seen in the Shannon water pouring theorem, where for increasing ambient 
noise and constant Transmission Loss, a continually greater SNR is required for signal 
reception.7 


7 Shannon, C.E., Weaver, W. (1949) A Mathematical Theory of Communication (Univ. of Illinois 
Press, Urbana, IL. 1949). 
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Figure 4. Environmental factor degrading channel Signal to Noise ratio for wind 
speeds ranging from 0 m/s to 25 m/s; variability due to wind noise is readily 

apparent. 
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III. EXPERIMENT 


Our experiment used a 40 node network deployed on the bottom of the oeean. To 
connect to the surface and send communications, a gateway buoy was moored at the 
southern end of the seafloor grid, as seen in Figure 5. 



Figure 5. The deployed network at COMEX; each circle represents one node, and 
the triangles represent gateway buoys. Each node is separated by 1-3 km (.5-2 

nm). 

This experiment was not purposefully designed to measure the performance of the 
underwater network with respect to wind. Therefore many measurements which would 
have helped to establish a firm relationship between range and wind speed such as 
measured ambient noise levels were not obtained. The gateway buoy recorded all 
network traffic that it could hear both migrating towards and away from it. That is, it 
recorded all message traffic, not just those messages that were addressed specifically to 
the gateway. Network performance was quantified by range to the furthest node each 
message could be heard traveling and decoded as it traveled away from or toward the 
gateway. The node number where the acoustic signal was last overheard was recorded, 

13 



















and used to find the range value using Table 1. Beeause of this, the ranges are diserete 
values, rather than having eontinuous range measurements. In short, the gateway buoy 
was the reeeiver, and all ranges were measured from there to the furthest transmitting 
node. 


Node 

Range in km 

Range in nm 

12 

0.76 

0.4 

13 

1.85 

1 

14 

2.6 

1.4 

15 

4.4 

2.4 

16 

5.6 

3 

17 

7.0 

3.8 

18 

8.5 

4.6 

19 

10.0 

5.4 

20 

11.1 

6 

21 

13.0 

7 

22 

13.7 

7.4 

23 

14.0 

8.2 

24 

15.5 

8.4 


Table 1. Ranges to eaeh node. 


A. MODEMS 

For this experiment. Benthos aeoustie modems were used in eaeh network node, 
as well as in the gateway buoy. To pass a message between nodes, the sending node first 
transmits a Request to Send (RTS) paeket. The reeeiving node responds with a Clear to 
Send (CTS) paeket. When that is reeeived, the sending node transmits the data of the 
message in eoneatenated 256 byte sub-paekets. If any of these paekets are eomipted, the 
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receiving node sends a Selective Repeat Request (SRQ) packet specifying which data 
sub-packet must be resent. 


Figure 6. 


Subsurface Float 
(24 pounds positive) 


(14 pounds negative) 


Telesonar Modem 


Acoustic Release 


Clump Weight 
(60 pounds negative) 



2m 


1m 


.5m 


1m 


.5m 


A schematic of one node, including the acoustic modem. 8 


B. SOURCE LEVEL AND FREQUENCIES 

The source level of each acoustic modem was set to 179 dB, and used 128 tones 
for l-of-4 MFSK. Each operated in the frequency band of 9-14 kHz. 

C. DATA 

The data were passed through the network in one of two types of packets: data 
and utility. 


Type of Packet 

Utility 

Data 

Raw Bit Rate 

2400 bps 

2400 bps 

Information Bit Rate 

150 bps 

800 bps 

Packet Size 

9 bytes 

Variable 


8 Rice, Joseph A. Seaweb Acoustic Communication and Navigation Networks from Proceedings of the 
International Conference “Underwater Acoustic Measurements: Technology and Results” Heraklion, Crete, 
Greece 28* June-U‘ July 2005. 
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Number of Packets 


656 


463 


Table 2. Comparison between Utility and Data Packets during experiment. 

1. Utility Packets 

Utility packets are transmitted at a bit rate of 150 b/s. Though both packet types 
have the same raw bit rate, the information bit rate of utility packets is slower because of 
additional coding. Theoretically, utility packets should be more reliable at longer ranges 
because of the larger amount of error correction and redundancy encoded within their 
data. Their size is quite small at 9 bytes. They are used for network maintenance 
messages, for instance. Clear to Send or Request to Send packets. The utility packets 
were used much more frequently during the experiment, providing us more data points. 

2. Data Packets 

Data Packets are used for actual data being sent through the network. These are 
larger in size at approximately 256 bytes, and are sent at 800 b/s. These typically 
contained 4 lines of text messaging. 
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IV. EXPERIMENT ENVIRONMENTAL EACTORS 


A. GEOMETRY 

The depth of the oeean at the gateway buoy was approximately 220 meters, 
gradually inereasing with range away from the buoy mooring. The nodes did not rest on 
the oeean bottom, but floated 3 meters above. Additionally, the modem transdueer on the 
gateway buoy was not at the surfaee; it hung 30-40 meters below. The distanee between 
the gateway node and the nodes of interest in the grid ranged from approximately 2 to 16 
km. 

To best approximate Transmission Loss, the unique bathymetry of the experiment 
should be modeled. Additionally, the loeal environment and sound speed profde should 
be represented in this model. 

B. SOUND SPEED PROFILE 

The sound speed profde of the oeean was not eonstant at the network, but varied 
throughout the eourse of the experiment as seen in Figure 7. The experimental sound 
speed profdes were taken using Conduetivity, Temperature, and Depth (CTD) sensors 
and Expendable Bathythermograph (XBT) measurement deviees, by members of the 
Naval Oeeanographie Offiee. 
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Figure 7. Combined SSP during the eourse of the experiment. 

A eharaeteristic sound speed profile from the experiment is included in Figure 8. 
Sound arcs toward the slower sound speed rather than traveling in a straight line. Notice 
the slight upward refraction in the mixed layer, extending to approximately 75 m. The 
mixed layer depth was not constant throughout the experiment. The main difference was 
the raising and lowering of the thermocline (the region of strong gradient just below the 
mixed layer). The change in the depth of the mixed layer was as much as 50 m. The 
region below the mixed layer was a downward refracting environment. 
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Figure 8. A simplified sound speed profile during the experiment, eompiled from 
approximately 10 separate SSPs, eourtesy of Paul Baxley, SPAWAR San Diego. 

C. TRANSMISSION LOSS VS. RANGE 

The best method of approximating Transmission Loss is in a physios based model 
whioh takes into oonsideration bathymetry, sound speed profiles, surfaoe roughness, and 
bottom type. There are several options available; however, a range-independent Kraken 
model is used. This is a normal mode propagation model. The assumptions are that the 
bottom is flat, the souroe depth is 217 m, the reoeiver depth is 50 m, the frequenoy is 11.5 
kHz (midway between 9 and 14 kHz), and the sound speed profile is that given in Figure 
8. The modeling was done by Paul Baxley of SSC San Diego. Five eases are modeled, 
eaeh with different assumptions about the boundary eonditions of the oeean ehannel. TL 
vs. range is plotted for eaeh ease in Figures 9-12. In eaeh figure, the dashed line is the 
plot of purely geometrie TL, due to a eombination of spherieal and eylindrieal spreading 
as shown in Figure 1. Figure 9 shows the TL ealeulated from the Kraken model without 
surfaoe or bottom losses. Note that at very short ranges the model prediots oonsiderably 
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higher transmission losses than the simplified geometrieal TL. At longer ranges the 
model also shows greater losses. At least part of this diserepaney ean be explained by the 
faet that the Kraken model ineludes losses due to absorption. 




Figure 9. 



Kraken Model of Transmission Loss vs. Range, with no surfaee or bottom 

losses. 
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Figure 10. Kraken model of Transmission Loss vs. Range for a lossy bottom and loss 

of 3 dB/bounee from surfaee refleetion. 
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Figure 11. Kraken model of TL vs. range for a lossy bottom and loss of 6 dB/bounee 

from surfaee refleetion. 








Figure 12. Kraken model of XL vs. range for a lossy bottom and loss of 9 dB/bounee 

from surface reflection. 

Figure 10 shows a difference of almost 10 dB greater TL when surface and bottom 
losses are introduced. As expected, Figure 11 shows even greater losses when the surface 
loss is 6 dB per bounce, which corresponds to four times the power loss for each bounce 
than for no surface losses at all. Finally, Figure 12 displays even worse conditions, with 
a lossy bottom and 9 dB loss per surface bounce. Using the Kraken model for different 
losses demonstrates the variability in TL for more complicated models rather than when 
using a simplified model such as geometric spreading. It also shows that the TL could be 
considerably higher than that calculated from simple geometrical spreading with 
absorption. 

D. RAY TRACING 

This was a downward refracting environment. Using sound speed profiles and 
local environment bathymetry, we can create the chart of possible ray paths from the 
seafloor transmitter (network nodes) into the ocean and eventually to the receiver 
(gateway buoy) displayed in Figure 13. 
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Figure 13. Ray traces for sound speed profde taken near the gateway buoy. The 
gateway buoy transducer depth was approximately 30-40 m, and bottom bounce 

rays are excluded. 
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V. EXPERIMENTAL WIND SPEEDS AND COMMUNCATION 

RANGES 


A. WIND VS. TIME 

Two sets of wind speed vs. time data are examined for this experiment. The first 
was forecast wind data, from COAMPS™ WESTPAC 27 km operational model run, 
provided by the Fleet Numerical Meteorology and Oceanography Center. This was the 
most complete data set, with a value for wind speed at every hour. Each forecast was 
based on an analysis performed at OOOOZ or 1200Z. For the purposes of this study, we 
were not able to determine if reports from ships in that area were incorporated in these 
analyses. However, these data do not agree well with the observed wind speed from the 
research vessel (RV, which was moored at the site of the network) during the experiment 
period. The other set of wind speed data is from the RV Bridge. The bridge data were 
recorded three times a day, each day. Figure 14 shows the wind data during the period of 
the experiment. 
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Wind Speed vs Day 
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Figure 14. Wind speed vs. time for all three measurements of wind data. The 

variability in wind speed over the eourse of the experiment provides the basis for 

this experiment. 


B. RANGE VS. TIME 

As stated previously, the measured eommunications ranges are diserete values, 
based on the distance to the furthest node heard during the sending of network packets. 
The gateway buoy recorded the nodes that were communicating, and the time at which it 
took place. The ranges corresponding to only the most distant overheard node are plotted 
against time to show network performance as the experiment continued. The 
communication from all closer nodes was ignored when plotting range, and all nodes are 
assumed to be communicating at an equal rate. However, since the undersea network 
was not in continuous use there are breaks in time where no traffic is observed. The 
following figure shows utility packet range vs. time. 
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utility Packet Range vs Tinne 



Figure 15. Utility Packet Range vs. Time 


Wind Speed vs. Day for RV data 



Figure 16. RV wind speed vs. time plot, for comparison to range vs. time 


However, there are ambiguous data in Figure 15. Since not all messages were 
destined for the farthest node, it is possible that some of the transmissions which are 
plotted as if they represent the farthest range attainable are actually terminating at the 
message destination, and not the limit of the gateway buoy acoustic range. If a message 
was addressed to a node at a range less than the maximum range the gateway could listen. 
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then it would stop before the limit of the gateway hearing was reaehed. That is, the data 
are incorrect because the maximum distance over which that signal was heard is shorter 
than it would have been if the packet had been destined for a more distant node. Culling 
those ambiguous data was not performed. The average range for utility packets was 6.77 
km, and the maximum range was 15.7 km. 


Data Packet Range vs Time 

^6r 

14 - 

12 - 


10 


4> 

OS 

c 


8 - 



2 - 


gl_I_I_I_I_ I _I_I_I_ I_» 

1 23456789 10 11 

Day 


Figure 17. The plot of Data packet range vs. time shows network performance. 
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Figure 18. RV wind speed vs. time for comparison to range vs. time 
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The data packet range and time plot is created in the same manner as the utility 
packet range vs. time plot, and has the same ambiguous data. For data packets, the 
average range was 6.33 km and the maximum range was 15.7 km, as well. 

C. RANGE VS. WIND SPEED 

Because there are far fewer wind speed vs. time data points than range vs. time, 
there is poor resolution in plotting wind and range against each other. Range was 
assigned a corresponding wind speed value for the wind speed observed prior to the time 
that each range was taken. This is then plotted with range as a function of wind speed, in 
Figure 19. 


Weighted Utility Packet Range vs. Wind Speed 



Figure 19. Weighted utility packet range vs. wind speed, with a 50% line using RV 

wind data. 


The plot of utility packet range and wind speed is weighted to show how often a 
particular range happened for a given wind speed. A larger circle denotes greater 
occurrence. The total number of all utility packets overheard at each wind speed was 
recorded. The 50% success line in Figures 19-22 represents the node for which half of all 
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utility packets were overheard either at or further than that node. For example, when the 
wind speed was 14m/s, half of the utility paekets overheard originated from the node 
about 2km away or from more distant nodes. Wind data for Figures 19 and 21 were 
observed at the RV Bridge, while the wind data for Figures 20 and 22 are from forecast 
data. Using the forecast wind data, a weighted oeeurrence plot of range and wind speed 
appears slightly different. The same processing is applied to the data packet ranges. 


Weighted Utility Packet Range vs. Forecast Wind Speed 
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Figure 20. Weighting the forecast wind speed vs. range shows a less organized plot, 

with less elear potential trends. 
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Weighted data packet range vs. wind speed for RV wind data. Larger 
circles are for more occurrences, and 50% success lines drawn. 

Weighted Data Packet Range vs. Forecast Wind Speed 
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Weighted data packet range vs. forecast wind speed; again, the trends 
between RV and forecast data differ. 
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VI. EXPERIMENTAL LINK BUDGET 


Using the information that has been experimentally eolleeted, we now revisit the 
link budget to ereate a new relationship between range and wind speed. Solving 
Equation 2.1 for TL yields: 

TL = PSL - NSL - SNR (6.1) 

The only term dependent on range is TL. Assuming NSL is dominated by wind noise, 
range is now determinable for any known wind speed. However, the relationship 
involves a transeendental equation, sinee: 

TL = 10 logio (r,) +10 logio (r) + ar (6.2) 

Substituting equation 6.2 into 6.1 gives: 

10 log 10 (r) + ar = PSL - NSL - SNR -10 loglo (r,) (6.3) 

Returning to ehapter 2 and using the relationship for NSL^md in equation 6.3 gives: 

10 logio (r) + ar = P5L - SNR - 10logio {r,) - {1.5yfw + C) (6.4) 

where C is a eonstant for a given frequeney. Thus we have a relationship for wind speed 
and range; however range eannot be found explieitly, save for at shorter ranges when 
absorption ean be ignored. (For signals at our eenter frequeney, 11.5kHz, the loss from 
absorption is ealeulated as 1.3dB/km.) Solving for the range in that ease yields: 

(-7.5Vw+rai-SA'R-101og,o(r,)-c) 

r = 10 ^ (6.5) 

Admittedly, this simplified model does not reflect the complete dependence of the 
communication range on wind speed since Transmission Loss will also dependent upon 
wind due to increased surface scattering losses. Over short ranges, however, the losses 
due to surface bounces are expected to be small in comparison to geometrical spreading 
and absorption. 

We can see from this simplified model that the expected communication range 
decreases with increasing wind speed. The limited data set available supports this 
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relationship. With a new and better data set, it may be possible to establish a more 
definitive relationship for network performanee prediction based upon wind speed. 
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VII. CONCLUSIONS 


This thesis intended to examine the eorrelation between eommunieations range 
and wind speed. To aehieve this, we ereated an example link budget as a benehmark for 
hypothetieal network performanee, and showed that theore tie ally, inereased wind speed 
degrades network performanee. Next, we examined experiment data and used range vs. 
wind speed plots to quantify eommunieations range against wind speed. Finally, we 
revisited the link budget, and used the experimental data to show the degradation of 
network performanee for inereasing wind speeds. Unfortunately, this experiment was not 
intended to test this hypothesis in a eontrolled and statistieally signifieant manner, and 
our result therefore eannot be eonsidered eonelusive. 

Although the data show a rough correlation between wind speed and 
communications range, a more exact relationship between the wind speed and range 
could not be established. A controlled experiment should be carried out to obtain more 
conclusive data. That experiment should have several improvements over the last. First, 
wind speed and sea conditions must be sampled with rigor. Repeated sound speed profde 
measurements should be taken to resolve any periodic changes of the mixed layer and 
thermocline. Making sure that many messages are sent within the network for each wind 
speed is a must. This will create many more data points for analysis. Ambient noise 
should be measured in both the mixed layer and in the thermocline. Finally, a directional 
receiver for measuring vertical arrival angles to verify the ray paths would be helpful in 
better understanding the multipath structure. Some things should be retained from the 
earlier experiment. The fact that the network stretched well beyond the hearing of the 
gateway buoy allowed for data to be recorded with no fear of acoustic signals 
propagating beyond the network while within hearing of the gateway buoy. Also, a time 
of year where significant wind speed changes are the norm is a must, to ensure a wide 
range of wind speeds. 
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Using a calibrated link budget model is a useful tool for network design. It allows 
for simple network deployment no matter what the eonditions or loeal environment. In 
order for the Department of Defense to exploit aeoustie eommunieations teehnology, 
more researeh is neeessary. 
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APPENDIX: MATLAB CODE 


A. MAIN LINK BUDGET PROGRAM 

%Link Budget Equation 
%Thesis Version 3 
elear all; ele; elose all; 

%initialize parameters 

[f, R, depth, ws, D] = initialize; 

% Transmission Loss (TL) 
dispCealeulating TL') 

[TL, TLa, TLs, alpha 1] = transmissionloss(R, f, depth); 

% Noise Level (NL) 
dispCealeulating NL') 

[NLwind, NLship, NLturb, NLtherm, NLtotal, NL] = noiselevel(R, f, ws, D); 

%Souree Level 

dispCealuelating souree level'); 

[SL] = soureelevel(R, f); 

%Link Margin 

dispCealeulating the link margin') 
linkmargin = SL + TL + NL; 

%Caleulate the 'Noise' to be overeome 
dispCealeulating the Noise to be overeome') 
relSNR = NL + TL; 
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%PLOT OUTPUTS 

dispCgenerating plots') 

% %TL plots 

allplots('TLplots', R, f, TL, TLa, TLs, alphal) 

% 


NL) 


% %NL plots 

allplots('NLplots', R, f, ws, D, NLwind, NLship, NLturb, NLtherm, NLtotal, 


% 


% % 




% %Relative SNR 

allplots('relSNRplots', R, f, ws, D, relSNR) 
% 

% He***** He**** **=(==[=* He* 


% %Link Margin plots 

allplots('linkmarginplots', R, f, linkmargin, depth) 

% 


B. INITIALIZE SUBROUTINE 

function [f,R,depth,ws,D,Pe,nu,L] = initialize 

%Initialize Data 

!lc sic !lc sic !lc ![C 5lc ![C !lc sic !lc sic !lc sic 5lc ![C !lc sic !lc sic !lc ![C 5lc * !lc sic 

clear all; clc; 

%Deline Range and Frequency Limits: 
%frequencies of interest (Hz) 
fupper = 100000; 
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f_lower = 0; 
f_res = 50; 

%ranges and depth (m) 
depth = 220; 
range = 1500; 
rangeres =10; 

%Generate Frequency and Range vectors: 

% BW = f_upper - f_lower; 

% 

% stop = ceil(loglO(f_upper+l)); 

% for i = l:stop; 

% f(:,i) = [ (l*10^(i-l)): (f_res*10^(i-l)): (9.99*10^(1-1)) ]’; 

% end 

% f = reshape(f,9*i/f_res,l); 

% begin = find(f==f_lower); 

% finish = find(f==f_upper); 

% 

% f = f(begin: finish); 

f=(f_lower:f_res:f_upper); 

f=f; 

R= [0:range_res:range]; R(l)=l; 

%Deline noise related parameters 

%wind speed in m/s 

ws=25; 

%shipping activity low (D=0), med (D=0.5), high (D=l) 

D=0.5; 

C. TRANSMISSION LOSS SUBROUTINE 

function [TL, TLattenuation, TLspreading, alphal] = transmissionloss(R, f, depth) 

% TRANSMISSION LOSS 
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%TL due to Spreading - Matrix (Range x Freq) 

Rindexl = find(R==depth); %max range of spherieal spreading 
%Rindex2 = lind(R==depth); %max range of transitional spreading 
RindexS = find(R==depth); %max range of eylindrieal spreading 
TLspreadingl = 20*logl0(R(l;Rindexl)); 

TLspreadingS = (10*logl0(R(Rindexl+l;end))-10 *loglO(R(Rindexl))) + 
TLspreadingl (end); 

TLspreadingtotal = [TLspreadingl TLspreadingS]; 

TLspreading = TLspreadingtotal; 
for findex = 2;length(f) 

TLspreading = [TLspreading; TLspreadingtotal]; 
end 

%TL due to Attenuation - Matrix (Freq x Range) 

%using eqn from Uriek p 108 

%note this is given in dB/kyd with f in kHz therefore need to eonvert 
fkhz = f/1000; 

alphal = ((0.1L*(fkhzA2))./ (1 + (fkhzA2))) + ((44 A (fkhz A2))./ (4100 + 
(fkhzA2))) + ((3.0A(10A(-4))).*(fkhzA2)) + 0.0033; 

for Rindex = l:length(R) 
for findex = 1 :length(f); 

TLattenuation( findex,Rindex) = alphal(findex).*R(Rindex).*le-3; 
end 
end 

%Total Transmission Loss Matrix 
TL = TLattenuation + TLspreading; 

TL=-TL; 
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D. NOISE LEVEL SUBROUTINE 


function [NLwind, NLship, NLturb, NLtherm, NLtotal, NL] = noiselevel(R, f, ws, 

% NOISE LEVEL 

fkhz = f./lOOO; 

%NL due to Surface Waves Option 2 frin 

NLwind = ( 50 + (7.5.*(ws.^0.5)) + (20.*logl0(fkhz)) - (40.*logl0(fkhz + 0.4)) ); 
%NL due to Shipping 

NLship = ( 40 + (20.*(D-0.5)) + (26.*logl0(fkhz)) - (60.Nogl0(fkhz + 0.03))); 


%NL due to Turbulence 
NLturb = ( 17 - (30.*logl0(fkhz)) ); 

%NL due to thermal noise in the receiver 
NLtherm = (-15 + (20.*logl0(fkhz)) ); 

%Rain 

%Industrial Noise 
%Biologics 

%Total Noise Level 

NLtotal = 10.*logl0( (10.^(NLwind./10)) + (10.^(NLship./10)) + 
(10.^(NLturb./10)) + (10.^(NLtherm./10))); 

NL = NLtotal; 

for Rindex = 2:length(R) 

NL = [NL NLtotal]; 
end 

NL=-NL; 
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E. SOURCE LEVEL SUBROUTINE 


function [PSL] = sourcelevel(R, f) 

% Source Level 

=1= He =1= He =1= =[c ** =lc He =1= He =1= =1= ** ****=[=* *=1= **=1=* *=1= * 

%SL = 171.5 + 10.*logl0(Pe) + 10.*logl0(nu); 

SL=179 

SL = SL.*ones(length(f),length(R)); 

PSL= SL-10*log(32); 

F. PLOTS SUBROUTINE 

funetion varargout = allplots(plots,varargin) 
switch plots 
case 'TLplots' 

% When the function allplots is called in linkbudgeteqnv4.m, this 
% section plots the Transmission Loss 

% It takes in the Range (R), Frequency (f), Total Transmission Loss (TL), 
% Transmission Loss due to attenuation (TLa), Transmission Loss due to 
% spreading (TLs) and the Attenuation coefficient 

R=varargin{l}; f=varargin{2}; TL=varargin{3}; TLa=varargin{4}; 
TLs=varargin {5 }; alpha 1 =varargin {6 }; 

%logarithmie plot Spreading Transmission Loss (in dB) vs Range (in m) 

figure(l) 

semilogx(R,TL(l,:)); title('Total Transmission Loss'); grid on; 
xlabel('Range (m)'); ylabel('dB'); 

%logarithmic plots the attenuation coeffieient vs freq (in kHz) 
figure(2) 

semilogy(f/1000,alphal); 

%title('attenuation eoeff); 
grid on; 

xlabel('Freq (kHz)'); ylabel('dB/Km'); 
xhm([l 100]); 
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%plots a 3-dimensional graph of Range (in m) vs freq (in kHz) vs 
%Transmission Loss due to attenuation (in dB) 

% figure(3) 

% poolor(R,f/1000,TLa); 

% shading interp; 

% eolorbar; 

% hold on; 

% [os, h] = oontour(R, f/1000, TLa, [-100:10:200], 'k'); 

% clabel(os); 

% xlabel('Range (m)'); ylabel('Freq (kHz)'); zlabel('dB'); 

% %title('TLattenuation in dB'); 

% grid on; 

% hold off; 

% grid on; 

%plots a 3-dimensional graph of Range (in m) vs freq (in kHz) vs Total 
%Transmission Loss (in dB) 

% figure(4) 

% poolor(R,f/1000,TL); 

% oaxis([-80 -40]) 

% shading interp; 

% eolorbar; 

% hold on; 

% [os, h] = oontour(R, f/1000, TL, [-200:20:200], 'k'); 

% %olabel(cs); 

% xlabel('Range (m)'); ylabel('Freq (kHz)'); zlabel('dB'); 

% title('TL total in dB'); 

% hold off; 

% grid on; 

case 'NLplots' 

%When allplots is called in linkbudgeteqnv4.m, this section plots the 
%Noise Level 

%It takes in the Range (R), freq (f), wind speed (ws). Directivity (DO, 
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%Noise Level due to wind (NLwind), Noise level due to shipping 
%(NLship), Noise level due to surfaee turbulenee (NLturb), Noise Level 
%due to thermal relaxation (NLtherm), Total Noise level (NLtotal), 

%Total Noise Level in matrix form (NL) 

R=varargm{l}; f=varargm{2}; ws=varargm{3}; D=varargm{4}; 
NLwind=varargin {5 }; NLship=varargin {6 }; 

NLturb=varargm{7}; NLtherm=varargm{8}; NLtotal=varargm{9}; 
NL=varargin {10}; 

%logarithmio plot of the various noise sourees and total noise (in dB) vs fireq 

(in kHz) 

figure(5) 

semilogx(f/1000,NLwind,'b',f/1000,NLship,'r',f/l 000,NLturb,'g',f/l 000,NLtherm,'e',f/100 
0,NLtotal,'k'); 

“/otitleCNL in dB'); 

grid on; 

xlabel('Freq (kHz)'); ylabel('dB'); 
legendCwind', 'shipping', 'turbulenee', 'thermal', 'total'); 
ybm([0 110]); 
xbm([0 100]); 

%gtext({sprmtf('wmd speed; %2.If m/s', ws),... 

%sprmtf('Shipping eoeff (D): %2.1f, D)}); 

% logarithmie plot of total Noise level (in dB) vs fireq (in kHz) 

% figure(6) 

% semilogx(f/1000,NLtotal); title('NLtotal in dB'); grid on; 

% xlabel('Freq (kHz)'); ylabel('dB'); 

% 


case 'relSNRplots' 

%when allplots is called from linkbudgeteqnvd.m this portion plots the 
%relative Signal to Noise Ratio 

% it takes in range (R), freq (f), wind speed (ws), directivity (D), 

% and relative Signal to noise Ratio (relSNR) 
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R=varargin{l}; f=varargin{2}; ws=varargin{3}; D=varargin{4}; 
relSNR=varargin {5 }; 

%plots the range in 1 km inerements vs relative SNR 
% figure(7); 

% Rlk=R(l 1); R2k=R(21); %R5k=R(51); R10k=R(101); R20k=R(201); 
% plot(f/1000,relSNR(:,l l),'k',f/1000,relSNR(:,21),'g'); 

% %,f,relSNR(:,51),’b’,f,relSNR(:,101),'e',f,relSNR(:,201),'r'); 

% grid on; 

% xlabel('Freq (kHz)'); ylabel('dB'); 

% %title('Total Envimmentals to Overeome'); 

% legend('l km', '2 km', '5 km', '10 km', '20 km'); 

% %gtext({sprintf('wind speed; %2.1f m/s', ws),... 

% % sprintf('Shipping eoeff (D): %2.1f, D)}); 

%plots a 3- dimesional graph of range (in m) vs freq (in kHz) vs 
%relative SNR (in dB) usin g a eolor range from -160 dB to -120 dB 
figure(8) 

poolor(R,f/1000,relSNR); 
shading interp; 
oaxis([-120 -80]); 
eolorbar; 
hold on; 

[os, h] = oontour(R, f/1000, relSNR,[-200;20;200],'k'); 
olabel(os); 

xlabel('Range (m)'); ylabel('Freq (kHz)'); zlabel('dB'); 
title('Total Environmental to Overeome in dB'); 
hold off; 

ease 'linkmarginplots' 

%when allplots is oalled in linkbudgeteqnv4.m, this part plots the link 
%budget margin 

%It takes in range (R), frequenoy (E), link margin (linkmargin), and 
%depth (depth) 
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R=varargin{l}; f=varargin{2}; linkmargin=varargin{3}; depth=varargin{4}; 


%this plots a 3-dimensional graph of range (in m) vs freq (in kHz) vs 

%linkmargin (in dB) 

figure(9) 

peolor(R,f/1000,linkmargin); 
shading interp; 
eolorbar; 
hold on; 

[os, h] = oontour(R,f/1000, linkmargin, [17.8 17.8], 'r'); 

%olabel(os); 

xlabel('Range (m)'); ylabel('Freq (kHz)'); zlabel('dB'); 
title('Link Margin in dB'); 
hold off; 
end 

G. PLOTTING WIND VS. TIME, RANGE VS. TIME, WIND VS. RANGE 

%u and V are oomponents of wind, and the time must be oonverted into 

%deoimal values for day; the text file oontains the foroast wind data 

load u2004110712.txt 

u=u2004110712(1:12,6); 

v=u2004110712(1:12,7); 

tautemp=u2004110712(1:12,5); 

tautemp=( 12-l-tautemp)/24-l-10; 
uvtemp=sqrt(u. ^2-l-v. ^2); 
load windseries 

wspd = [wspd' uvtemp(l: 11)']'; 
jdate = [jdate' tautemp(l:l 1)']'; 
save windseries wspd jdate 

%Prof Rioe's wind observations 
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ricewind=[15.4 15.4 10.3 15.4 15.4 10.3 7.7 5.1 7.7]; 
ricetime=[l;9]; 


load Narra_wind.mat 
loadNarra time.mat 


figure(2) 

p lot(ricetime ,ricewind, 'r:') 

titleCWind Speed vs Day'), xlabel('Day'), ylabel('Wmd Speed (m/s)') 
hold on 

plot(jdate,wspd,'-.') 
p lo t(Narra_time,Narra_wmd,'g—') 

%need to read in the range data from exeel 

%need file 'thesis data.xls' and sheets 'utility range' and 'data range' 

Utility=xlsread('Thesis Data.xls','Utility Paekets'); 
Data=xlsread('Thesis Data.xls','Data Paekets'); 

load dat.mat 
load util.mat 

figure(lOl) 
subplot(3,l,l) 
plotfjdate, wspd) 

title('Foreoast Wind Speed Data vs. Time') 
ylabel('Wind Speed (m/s)') 
subplot(3,l,2) 
plot(rieetime,rieewind) 

title('Prof Riee Observed Wind Speed Data vs. Time') 

ylabel('Wind Speed (m/s)') 

subplot(3,l,3) 
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plot(Narra_time,Narra_wind*0.514) 
title('Narragansett Observed Wind Speed Data vs. Time') 
ylabel('Wind Speed (m/s)'),xlabel('Time (days)') 

figure(3) 

seatter(util(:, 1) ,util(; ,2) ,4) 

xlabel('Day'),ylabel('Range (km)'),title('Utility Paeket Range vs Time') 
figure(4) 

seatter(dat(:, 1 ),dat(: ,2),4) 

xlabel('Day'),ylabel('Range (km)'),title('Data Paeket Range vs Time') 

%Need to get both wind speed and range on the same time seale for foreeast 
%data 

%double for loop to synehronize time seales 
for i = 1:463 
forj = 1:220 

if jdate(j) >= dat(i,l); 
windvalue(i) = wspd(j-l); 
break 
end 
end 
end 

%repeat loop for utility paekets 
for i = 1:656 
forj = 1:220 

if jdate(j) >= util(i,l); 

windvalue2(i) = wspd(j-l); 
break 
end 
end 


end 
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save windvalue 
save windvalue2 

%get range and windspeed on same time seale for Prof Riee's observations 
for i = 1:656 
for j = 1:9 

if ricetime(j) >= util(i,l); 
windvalue3(i) = ricewind(j-l); 
break 
end 
end 
end 

for i = 1:463 
for j = 1:9 

if ricetime(j) >= dat(i,l); 
windvalue4(i) = ricewind(j-l); 
break 
end 
end 
end 

%Need to get both wind speed and range on the same time seale for M/V 
%Narragansett observations 
%double for loop to synehronize time scales 
for i = 1:463 
for j = 1:30 

if Narra_time(j) >= dat(i,l); 

windvalueN_D(i) = Narra_wind(j-1); 
break 
end 
end 
end 
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%repeat loop for utility packets 
for i = 1:656 


for j = 1:30 

if Narra_time(j) >= util(i,l); 
windvalueN_U(i) = Narra_wind(j-1); 
break 
end 
end 
end 

save windvalueN_U 
save windvalueN D 


figure(5) 

seatter(windvalue,dat(:,2),5) 

title('Data Paeket Range vs Foreeast Wind Speed') 

ylabel('Data Paeket Range (km)'),xlabel('Wind Speed (m/s)') 

figure(6) 

seatter(windvalue2 ,util(: ,2), 5 ) 

title('Utility Paeket Range vs Foreeast Wind Speed') 

ylabel('Utility Paeket Range (km)'),xlabel('Wind Speed (m/s)') 

figure(7) 

p lotyy (dat(:, 1) ,dat(: ,2), dat(:, 1) ,windvalue) 
figure(8) 

plot(rieetime,rieewmd,util(:,l),util(:,2),'*') 

title('Utility Paeket Range and Observed Wind vs Time') 

xlabel('Day') 


figure(9) 
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p lot(ricetime ,ricewind, dat(:, 1) ,dat(: ,2),' * 

title('Data Packet Range and Observed Wind vs Time') 

xlabel('Day') 

figure(lO) 

seatter(windvalueN_D,dat(:,2),'*') 

title('Data Paeket Range vs. M/V Narragansett Wind Speed') 
ylabel('Data Paeket Range (km)'),xlabel('Wind Speed (m/s)') 

figure(ll) 

semilogy(windvalueN_U,util(: ,2),' *') 

title('Utility Paeket Range vs. M/V Narragansett Wind Speed') 
ylabel('Utility Paeket Range (km)'),xlabel('Wind Speed (m/s)') 

H. WEIGHTING WIND VS. RANGE AND PLOTTING 

%program to weight range vs. wind speed plot 

elear all; ele; elose all; 

load util.mat 
load dat.mat 
load windvalueN_U 
load windvalueN_D 
load Narra_time 
load Narra_wind 

weightU=zeros(13,9); 

weightD=zeros(13,9); 

ranges=[2 5 7 12 15 19 23 27 30 35 37 41 42]; 

windspeed=[Narra_wind(12) Narra_wind(8) Narra_wind(16) Narra_wind(18) 
Narra_wind(9) Narra_wind(15) Narra_wind(l) Narra_wind(4) Narra_wind(2)]; 
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for a=l:13 
for b=l;9 
for c=l:656 

if ranges(a)==util(c,2) 

if windspeed(b)==windvalueN_U(c) 
weightU(a,b)=weightU(a,b)+l; 
end 
end 
end 
end 
end 

figure(12) 

axis([2 16 0 45]),title('Weighted Utility Paeket Range vs. Wind Speed') 
ylabel('Utility Paeket Range (nm)'),xlabel('Wind Speed (m/s)') 
for a=l;13 
for b=l:9 
if weightU(a,b)>0 
hold on 

soatter(windspeed(b),ranges(a),weightU(a,b)*4,'b','fdled') 

end 

end 

end 

seatter(windvalueN_U,util(: ,2), 1) 

ninetylineU=[ranges(9) ranges(8) ranges(8) ranges(8) ranges(8) ranges(6) 
ranges(4)]; 

fiftylineU=[ranges(7) ranges(6) ranges(6) ranges(6) ranges(4) ranges(4) 
ranges(2)]; 

windspeed2=[windspeed(3) windspeed(4) windspeed(5) windspeed(6) 
windspeed(7) windspeed(8) windspeed(9)]; 

p lot(windsp eed2 ,ninety lineU, 'r') 

plot(windspeed2,fiftylineU,'g') 
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for a=l:13 
for b=l;9 
for c=l:463 

if ranges(a) ==dat(c,2) 

if windspeed(b)==windvalueN_D(c) 
weightD(a,b)=weightD(a,b)+1; 
end 
end 
end 
end 
end 

figure(13) 

axis([2 16 0 45]),title('Weighted Data Paeket Range vs. Wind Speed') 
ylabel('Data Paeket Range (nm)'),xlabel('Wind Speed (m/s)') 
for a=l;13 
for b=l:9 
if weightD(a,b)>0 
hold on 

soatter(windspeed(b),ranges(a),weightD(a,b)*4,'b','fdled') 

end 

end 

end 

seatter(windvalueN_D,dat(:,2),l) 

ninetylineD=[ranges(9) ranges(6) ranges(8) ranges(8) ranges(8) ranges(5) 
ranges(5)]; 

liftylineD=[ranges(6) ranges(6) ranges(6) ranges(4) ranges(4) ranges(4) 
ranges(2)]; 

p lot(windsp eed2 ,ninety lineD, 'r') 
plot(windspeed2,fiftylineD,'g') 

%repeat weighting, for foreeast data 
load windseries 
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weightUF=zeros( 13,220); 
weightDF=zeros( 13,220); 

wspd=sort(wspd); 

for a=l;13 
for b= 1:220 
for c=l:656 

if ranges(a)==util(c,2) 

if wspd(b)==windvalue2(c) 
weightUF (a,b)=weightUF (a,b)+1; 
end 
end 
end 
end 
end 

for a=l:13 
for b= 1:220 
for 0=1:463 

if ranges(a)==dat(c,2) 

if wspd(b)==windvalue(o) 

weightDF (a,b)=weightDF (a,b)+1; 
end 
end 
end 
end 
end 

ninetylineUF=[ranges(9) ranges(lO) ranges(6) ranges(8) ranges(9)]; 
fiftylineUF=[ranges(6) ranges(8) ranges(4) ranges(5) ranges(6)]; 
lines=[3 4 5 6 7]; 
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figure(14) 

axis([2 9 0 45]),title('Weighted Utility Packet Range vs. Foreeast Wind Speed') 
ylabel('Utility Paeket Range (nm)'),xlabel('Wind Speed (m/s)') 
for a=l;13 
for b= 1:220 

if weightUF(a,b)>0 
hold on 

soatter(wspd(b),ranges(a),weightUF(a,b)*4,'b','fdled') 

end 

end 

end 

plot(lines,ninetylineUF,'r') 

plot(lines,fiftylineUF,'g') 

ninetylineDF=[ranges(9) ranges(lO) ranges(6) ranges(8) ranges(8)]; 
fiftylineDF=[ranges(6) ranges(8) ranges(4) ranges(4) ranges(6)]; 

figure(15) 

axis([2 9 0 45]),title('Weighted Data Paeket Range vs. Foreeast Wind Speed') 
ylabel('Data Paeket Range (nm)'),xlabel('Wind Speed (m/s)') 
for a=l:13 
for b= 1:220 

if weightDF(a,b)>0 
hold on 

soatter(wspd(b),ranges(a),weightDF(a,b)*4,'b','fdled') 

end 

end 

end 

plot(lines,ninetylineDF,'r') 

plot(lines,fiftylineDF,'g') 
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